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Army  Needs  in  Reliability,  Maintenance 

and  Logistics 


>  Reduce  operations  and  maintenance  costs 

>  Increase  effectiveness  of  fleet  logistics 

>  Control  lifecycle  cost  and  also  use  it  in  design  and 
procurement 

>  Improve  availability;  schedule  maintenance 

>  Use  both  analytical  and  experimental  /  field  data  to 
estimate  reliability 
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Background 


Random  Process 
(Field) 


Time  t 
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Input 

Random 

Process 


Vehicle 


Output 

Random 

Process 


Random  Process  leads  to  Time-Dependent  Reliability 
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Research  Statement 
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Develop  methodologies  to  obtain  a  preventive 
maintenance  schedule  and  to  assess  and  improve  the 
reliability  /  durability  of  vehicle  systems  using 


Experimental  (field) 
“Expert”  opinion 


Previously  and 
currently  at 
TARDEC 


Predictive  tools  (physics-of-failure  data) 


Current  research 


5/24/2011 


UNCLASSIFIED 


4 


Overview 
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Part  1: 

Optimal  preventive  maintenance  schedule 

using  time-dependent  reliability  and  lifecycle  cost 

Part  2: 

Accelerated  testing  method  based  on 
importance  sampling  using  few  tests  which  run  for 
only  a  short  time 
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Part  1 :  Optimal  Preventive 
Maintenance  Schedule 
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What  is  Reliability? 
Cumulative  Probability  of  Failure 


m 
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Reliability  at  time  t  is  the  probability  that  the  system 

has  not  failed  before  time  t. 


such  that  g(x(t),t)<  0) 


Cumulative 
Prob.  of  Failure 


Instantaneous  Prob.  of  Failure 


Calculation  Methods  for  Fj{t) 

■  Maximum  Response  Method 

■  Niching  GA  &  Lazy  Learning  Local  Metamodeling 


■  Analytical 


■  MCS  /  Importance  sampling 

Fj  (t)  =  1  -  exp  [- j  A(t  )dt] 
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Definition  of  Lifecycle  Cost 
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Lifecycle  Cost  =  Production  Cost 


Quality 


ected  Variable 


[ 


Time-Dependent  System  Reliability 
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Definition  of  Lifecycle  Cost 
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C,(d,X,tf,r)=Cp{d,X)+  C,  (d,  X,  t„  )+  C,f  (d,  X,  f ,  r) 


/ 


Lifecycle 

Cost 


/ 


\ 


Production 

Cost 


Inspection 

Cost 


/ 

Expected 
Variable  Cost 


Final  time 


Interest  rate 

C*{d,X,tf,r)=  ]cP(t)e-r,f;(t)dt 


0 


I 


\ 


Cost  of  failure 
at  time  t 


PDF  of  time  to 
failure  time 


f;(0=p(3(e[o,(J  s  uch  that  g  (x(^ ),  t )  <  0) 
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Preventive  Maintenance  Schedule 
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Estimation  of  Time  for  Preventive  Maintenance 


max  t 


M 


d,Px  ’ax  Jm 

s.  t.  CL(d,  ’  ^)  — 

Fj  (d, X,tM)<l-R' (tM 

dL<d<d  u 
PxL  —  M^x  — 
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<  Ov  <  o 
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\ 

Time  for 
Maintenance 
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A  Roller  Clutch  Example 
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Random  Variables:  D,  d,  A 

Due  to  degradation: 

D->D(l-fcf) 

d  ^  d(l  -  kt) 

A  ->  A(l  +  kt) 

with:  k  =  2.5E-0A  mm/ year 
5/24/2011 


Constraints: 

Contact  angle  a  =  0. 1 1  ±  0.06  rad 

Torque  r>=3000  Nm 
Hoop  stress  crh  <  400  MPa 


gx{D,  d,  A)  =  0.05  —  cos  1 


(  D-d 


g2(D,d,A )  =  cos 


V  A  —  d 
/  D-d'' 


\ 


J 


<0 


V 


A  —  d 


—  0.17  <  0 


J 


g3(D,d,A)=3000-NL 


g4(D,d,A}  = 


Vv 


V  C1  J 


—  4E¥<o 

4  (D  +  d) 


N 

r-(.Y 

Y  Dd  'l 

s 

(b2  +  a2' 

2  71 

vciy 

K(D  +  d); 

A 

[b2-a2) 

-  400E06  <  0 
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Roller  Clutch:  Lifecycle  Cost 


where: 


CP  = 


7  r  0.75) 

f 

0.65) 

7  0.88^ 

3.5+ 

+ 

3.0  + 

+ 

0.5+ 

V  3(J/)  j 

V 

3(U  J 

V  3<ta  y 

7 


C,  =  20  F’.  (X,  t0 ) 

Scrap  cost/unit 


-fiifYt 


7 


=  10 
r  = 


Failure  cost/unit  (warranty  cost) 
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years 

3% 
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Roller  Clutch:  Reliability  vs  Time-to- 

Maintenance 

System  Reliability  vs  Time 

- Initial  design 

- Optimal  design 


0.5  - 


Maximum  t. 

limited  by 
side  consti'aints 


t  =  7.01  years 


t  =12.77  yeai's 


16  18  20 

Time  (years) 
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Roller  Clutch:  Pareto  Optimality 
between  Time-to-Maintenance  and  Cost 


Target  Lifecycle  Cost  vs  Time  for  Maintenance 
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Roller  Clutch:  Pareto  Optimality 
between  Time-to-Maintenance  and  Cost 
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Design  Variables: 

Hx  =  9  Pd  ’  Ma  }  °X  =  {°D  ’  ’  °A  } 


Side  Constraints: 

55.0973<  jUD  <55.4973  22.66  <  juD  <  23.06  101.49<  //A  <101.89 

0.04  <  cr p  <  0.08  0.03  <^<0.1  0.07  <  <rA  <  0. 1 1 3 


c[ 

18 

19 

20 
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25 

26 

27 

28 

Md 

55.4946 

55.4973 

55.4973 

553822 

55.4973 
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55.4973 

55.4973 
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Md 
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22.7735 
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22.8535 
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22.8208 
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22.8315 

22.8316 

Ma 

101.49 

101.49 

101.49 

101.49 

101.49 

101.49 

101.49 

101.49 

101.49 

101.49 

101.49 

0.08 

0.08 

0.0771 

0.0693 

0.0661 

0.0593 

0.054 

0.0496 

0.0423 

0.04 

0.04 

0.0639 

0.0543 

0.0481 

0.0449 

0.0407 

0.0368 

0.0334 

0.0306 

0.03 

0.03 

0.03 

°A 

0.1107 

0.0946 

0.084 

0.0763 

0.0701 

0.07 

0.07 

0.07 

0.07 

0.07 

0.07 
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Part  2:  Accelerated  Testing 
using  Importance 
Sampling 


5/24/2011 


UNCLASSIFIED 


Problem  Description 


Oakland 


LMVrKgTYI 


Random 

Variables 


■3  - 1 - 1 - 1 - 1 - 

0  200  400  600  800  1000 

Longitudinal  Distance,  ft 


Terrain 


Vertical  Accel.  (G) 

■ 

3 

R  2 
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d  i 
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<  -1 
•^3 

>  -2 

20  40  60  80  100 

Time 


Vehicle  speed  :  20  mph;  Mission  distance  :  100  miles 

Simulation  can  be  practically  performed  for  a 

short-duration  time 
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Our  Approach 


A  novel  MC-based  method  to  calculate  the  time- 
dependent  reliability  (cumulative  probability  of 
failure)  based  on : 

>  short-duration  data  and  an  exponential 
extrapolation  using  MCS  or  Importance 
Sampling  ( infant  Mortality) 

>  Poisson’s  assumption  ( Useful  Life) 
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Efficient  MC  Simulation  Approach 


Oakland 


i 


LM'.'rREnv 


5/24/2011 


UNCLASSIFIED 


19 


Poisson  Assumption 


Oakland 
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Fr(t^,,t)  =  l-(1-F' (t.  ))<T”' 


T  V  min 


nun 


where  : 


t 

m,  =  £[A7(/min  ,t)]=  J v+{t)dt  =  vm (t  -  tmm  ) 

\  in 

Number  of  out-crossings 


1 


(?)  =  lim 

Ar^0,Ar>0 


P[g(d, X,t)>On  g( d, X, t  +  At) <  o] 


At 


Out-crossing  rate 
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Quarter- Car  Model  on  Stochastic  Terrain 
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Constant  design  parameters: 
ms  =1000  kg 
mu  =100  kg 

Vehicle  speed  =  20  mph 

Random  Input  variables 

Damping,  bs  ~  N(7000, 14002) 
Stiffness,  ks  ~  J/V(40xl03,(4xl03)2) 


Random  Output  Process 

(Vertical  Acceleration,  G’) 


Random  Input  Process:  Experimental  Stochastic 
Terrain  from  Yuma  Proving  Grounds. 


Threshold  =  2G 
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Quarter- Car  Model:  Road 
Input  Random  Process  Characterization 
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AR(3)  model  was  identified  based  on: 


Autocorrelation 

Function 


Sample  Autocorrelation  Function  (ACF) 


T 

-?tTT 

t  1 1  It. 

Mi 

l 

IMvr 

°-40  10  20  30  40  50 

Lag 


Autocorrelation  of 


Residual  process 

Sample  Autocorrelation  Function  (ACF) 


• 

•  T  t  T, 

•  -  t 

•  •  T  ? 

Tt  .Ttt  f. 

*  1  1 

Mill  i*i* 

1  ill-*-! 

ii**  * 

'  4 

■0.2 - - - - - 

0  10  20  30  40  50 

Lag 


U-  =  1.2456 

v 


wM -0.2976  w._2  - 0.1954 


m._3  +  (0,  0.5 1322) 


Statistical  tests  were  performed  to  verify  the  model 
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Failure  Rate 


Quarter-Car  Model:  Results 

(Failure  Rate  Estimation  for  Threshold  =  2G) 
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0.15 
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0.05 
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f  '  Y 


|l  nV 


' 


Estimated  parameters: 


A0  =0.1708 
b  —  0.08 1 8 


0.5 


1 

Time,  sec 


1.5 


-  True  value,  MCS 

- Efficient  MCS 

Eitoonential 


extrapolation 

Poisson's  Assumption 


40  60 

Time 


Estimation  requires  short 
duration  MCS 


Exponential 

extrapolation 


X(t)  «  /l0e  bt 
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Quarter- Car  Model:  Results 

Cumulative  Probability  of  Failure  for  Threshold  =  2G 
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Efficient  MCS  (blue)  approach  is  close  to 
true  MCS  results  (red) 
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Principle  of  Importance  Sampling: 
Random  Variable  Case 
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4^^  Importance  Sampling  for  Random  Process 
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Instantaneous  Conditional  Probability  of  Failure: 

Pf  ( t{ )  =  1 6>(x ;  tt  )/x  (x ;  tt  )dx 


Q 


x  =  {x1,x2,...,xi}  where  xj  is  a  realization  of  R.V.  Xt  =  X  (tj ) 
pj  & )  =  f  <9(x ;  t. )  fxs  (x;  f,  )nfx 


Sampling  Joint  PDF 


^/fc) 

Z  ^x;r,  )4x,0  Z^(x,0 

/#,)  =  — 
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4^^  Importance  Sampling  for  Random  Process 

ns 

Z  Z®(x’0 

P/  (*,-)=— 
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A'sfe-,) 


7l(t, )  =  lim  — —  =  lim  "  1 


NfM 

X  ^(x>  ) 


A;->0  At  Ar^O  A?  •  /Vv  (t;1 ) 


C0(X,t:)  = 


fxiX’tj) 

fxs (x;h) 


:  Likelihood  ratio  at  t. 


Ns  (t,  | ) :  Safe  sample  points  at  t,_, 


N  fitf)  :  Number  of  failures  in  At  =  t,  -  t._ , 
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4^^  Importance  Sampling  for  Random  Process 
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Likelihood  ratio: 

oj(x 


Decorrelation  length  :  Maximum  number  of  lags  over 
which  realizations  of  xt  are  significantly  correlated 


*i  -  V  =  01  (Xi-1  -  M)  +  <t>2  (Xi-2  -  M)  + . +  <Pn  ) 


To  generate  sampling  PDF 


fx  (x)  = 


1 


(2n)kl1  IE,1/2 


exp 


lv-1 


(l  ±  _± 

--(x-p)  jl  (x-p) 

v  2 


A 


) 
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Importance  Sampling  for  Random  Process 
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Estimation  of  Safe  Sample  Functions 


Nf{tt) 

A.  (t)  =  lim  — ^ - 

a^oA  t-Nsit^) 


“Inflated”  response 
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Quarter- Car  Example 
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ut  =1.2456 


—  0.2976m(_2  —  0.1954m._3 


+  £;(0,0.51322) 


Original  PDF  =  0.51  Sampling  PDF  <r5  =  0.7 


0  20  40  60  80  100  0  20  40  60  80  100 

Time,  sec  Time,  sec 


The  sampling  PDF  results  in  more  failures 
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Sample  Autocorrelation 


Quarter- Car  Example 


Sample  Autocorrelation  Function  (ACF) 
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Failure  Rate,  A.(t) 


Quarter- Car  Example 


1.5 


x  10 


-3 


Stationary  Case 


l 
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MCS  1,500,000  sample^ 
IS  10,000  samples 
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X  10 


-4 


Quarter- Car  Example 
Stationary  Case 
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Failure  Rate 


Quarter- Car  Example 
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Non-Stationary  Case 


Threshold  =  2  g  Threshold  =  2.65  g 


0.014 

0.012 

0.01 

|  0.008 

2a 

|  0.006 
Hh 

0.004 

0.002 

0 


Time,  sec 
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Observations  /  Practical  Issues 


>  Analytical  methods  can  be  used  under  the  Poisson’s 
assumption 


>  IS  at  initial  time  may  need  a  few  thousand  output  sample 
functions 
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Ongoing  Work  Plan 
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>  Improve  the  current  accelerated  testing  method  based  on 
importance  sampling  so  that  only  5-10  tests  are  needed  (Q3) 

S  Characterize  the  “inflated”  output  random  process  in 
importance  sampling  using  “generalized”  Kriging  and 
MLE  and/or  time  series 

>  Demonstrate  the  accelerated  testing  methodology  using  the 
N-post  (or  4-post)  Reconfigurable  Road  Simulator  of  the 
Physical  Simulation  Laboratory  at  TARDEC  (Q3  and  Q4) 


5/24/2011 


UNCLASSIFIED 


36 


5/24/2011 


UNCLASSIFIED 


37 


Thanks  for  your 
attention  ! 

Q&A 

**Disclaimer:  Reference  herein  to  any  specific  commercial  company,  product,  process,  or  service  by  trade  name,  trademark, 
manufacturer,  or  otherwise,  does  not  necessarily  constitute  or  imply  its  endorsement,  recommendation,  or  favoring  by  the 
United  States  Government  or  the  Department  of  the  Army  (DoA).  The  opinions  of  the  authors  expressed  herein  do  not 
necessarily  state  or  reflect  those  of  the  United  States  Government  or  the  DoA,  and  shall  not  be  used  for  advertising  or  product 
endorsement  purposes.** 
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